质子辐照硼硅酸盐玻璃盖片的物理效应分析

作为航天器电源系统的重要组成部分,太阳电池需要更高的转换效率和可靠性以及更长的使用寿命。通过在太阳电池表面覆盖抗辐照玻璃盖片,可以增强太阳电池对粒子辐射的防护,延长太阳电池的服役寿命,使航天器获得可靠的能源供应。硼硅酸盐玻璃就是一种理想的太阳电池玻璃盖片材料。采用蒙特卡罗方法,结合SRIM软件模拟研究质子辐照硼硅酸盐玻璃的损伤物理机理。基于粒子与物质相互作用的理论以及基本公式,通过分析不同入射能量的质子在硼硅酸盐玻璃中的阻止本领、电离能损、位移能损、空位的产生情况,对辐照损伤的物理机制进行研究。结果表明:能量为30~120 keV的质子辐照损伤主要发生在硼硅酸盐玻璃表面;质子沉积、空位分布等均为Bragg峰型分布;电离能损是能量损失的主要部分,随入射能量的增加而增大,导致电子的电离和激发;位移能损在玻璃内部随能量降低而增大,导致硼、氧和硅等空位缺陷的产生;电离效应和缺陷的产生是硼硅酸盐玻璃色心形成的重要原因。     

Structural,electronic, magnetic properties and critical behavior of the equiatomic quaternary Heusler alloy CoFeTiSn

In this article, we study the exchange coupling interactions of the equiatomic quaternary Heusler alloy CoFeTiSn, using the two methods: Monte Carlo simulations and the ab-initio method. In a first step, we use the ab-initio calculations to investigate the structural, the electronic and the magnetic properties of this alloy under the GGA method. The analysis of the energy dependence on the lattice parameter a (Å) of the equiatomic quaternary Heusler alloy CoFeTiSn, is discussed for different atomic configurations. The ferromagnetic configuration is found to be the more stable one, with an optimal lattice parameter value 6.00 Å. On the other hand, the electronic structure results show that the compound CoFeTiSn exhibits a half-metallic character and a spin polarization of 100% at the Fermi-level. The total magnetic moment of this alloy is found to be equal to 2.00 μ_B which follows the Slater Pauling rule. Our results support the half-metallic behavior of the studied material. In order to complete this study, we reported the dependence of the critical transition temperature as a function of the parameter α of the equiatomic quaternary Heusler alloy CoFeTiSn. We showed that the critical temperature increases almost linearly with an increase of the values of the parameter α.     

241Am-Be中子参考辐射的蒙特卡罗模拟研究

研究用于校准场所中子剂量监测仪表的241Am-Be中子参考辐射场计量特性。采用蒙特卡罗方法模拟了空气自由中子参考辐射(FRNR),GB/T 14055规定的最小尺寸中子参考辐射(SRNR)和实际中子参考辐射(ARNR)中不同检验点处中子周围剂量当量率、散射中子占比和能谱分布特征。研究结果表明,空气对FRNR中的剂量率和能谱分布影响小,近似为理想中子参考辐射;采用5%含硼聚乙烯作屏蔽的最小尺寸SRNR可减少热中子,降低散射中子占比,影锥法不适用于小尺寸中子参考辐射中对散射中子的修正;ARNR中的散射中子更少、占比更低,影锥法所得散射中子占比与理论值基本一致。     

均匀磁场下碳离子笔形束的剂量变化分析及位置修正方法

分析外加均匀磁场对于碳离子笔形束剂量分布的影响,并考虑修正这种影响,为磁共振成像引导碳离子放射治疗的临床应用提供指导。本文利用蒙特卡罗方法模拟计算了不同能量碳离子笔形束在不同强度磁场下的剂量分布情况,发现垂直于碳离子束入射方向的均匀磁场对于碳离子笔形束射程缩短的影响很小,磁场对碳离子束的主要影响是引起束流横向偏转,特别是碳离子束布拉格峰位置的横向侧移。横向侧移程度与碳离子束的能量和磁场强度相关,根据模拟结果,得到了一个计算碳离子束布拉格峰在磁场中相对横向偏转的方程,并提出一种校正外加磁场引起的碳离子束布拉格峰横移的角度修正方法。这些结果可用于指导磁共振图像引导碳离子放射治疗计划系统的研发。     

Monte Carlo study of the critical and the compensation behaviors of the double perovskite Sr2CrIrO6

The phase diagrams and magnetic properties of double perovskite $S{r}_{2}CrIr{O}_6$ have been studied by using Monte Carlo simulation based on the heat bath algorithm. The ground-state diagrams of the compound $S{r}_{2}CrIr{O}_6$ have been calculated for different combinations of system parameters. The diagrams obtained are very rich and they give an idea of all the most stable configurations. The effects of the exchange interactions and the crystal field on the phase diagrams and magnetic properties of the system have been examined. A number of interesting phenomena have been observed such as the compensation temperature, the first and second order phase transitions, the critical triple point and the terminal critical point.     

基于蒙特卡罗方法的中国中长期煤炭需求预测

为避免传统预测方法的参数取值主观性问题,采用参数随机产生的蒙特卡罗方法预测中国中长期煤炭需求。首先分析了经济增长、能源结构和产业结构三个主要煤炭需求影响因素,并基于1980~2015年间各影响因素及煤炭消费的历史数据和最小二乘法的多元线性回归拟合煤炭需求方程。在此基础上,构建各影响因素的概率分布,采用蒙特卡罗方法模拟1981~2015年的煤炭需求,发现仿真结果可以较好拟合现实,可作为仿真预测的有效工具。结合经济新常态和能源结构调整的现状,控制参数取值范围进行蒙特卡罗仿真预测,结果显示,2016~2025年的煤炭需求呈先上升后下降趋势,并于2020年达到需求的峰值40.25亿吨,这些结果对于煤炭产业的科学决策有重要的作用。     

Statistical Physics Approach to Liquid Crystals: Dynamics of Mobile Potts Model Leading to Smectic Phase,Phase Transition by Wang–Landau Method

We study the nature of the smectic–isotropic phase transition using a mobile 6-state Potts model. Each Potts state represents a molecular orientation. We show that with the choice of an appropriate microscopic Hamiltonian describing the interaction between individual molecules modeled by a mobile 6-state Potts spins, we observe the smectic phase dynamically formed when we cool the molecules from the isotropic phase to low temperatures (T). In order to elucidate the order of the transition and the low-T properties, we use the high-performance Wang–Landau flat energy-histogram technique. We show that the smectic phase goes to the liquid (isotropic) phase by melting/evaporating layer by layer starting from the film surface with increasing T. At a higher T, the whole remaining layers become orientationally disordered. The melting of each layer is characterized by a peak of the specific heat. Such a succession of partial transitions cannot be seen by the Metropolis algorithm. The successive layer meltings/evaporations at low T are found to have a first-order character by examining the energy histogram. These results are in agreement with experiments performed on some smectic liquid crystals.     

Kinetic study of hydrogen sulfide decomposition on Pt(111) surface

In this work, kinetic of H₂S conversion to H₂ molecule on the surface of Pt(111) is studied using kinetic Monte Carlo simulation. The results of simulation were fitted to the experimental temperature-programed desorption spectra. The good agreement between the empirical and the simulated data confirms the proposed mechanism and kinetic data (activated energies and pre-exponential factors). The influence of variables such as temperature and concentrations of H₂S and H₂ on the overall results of hydrogen production is studied. The condition is proposed in which the best yield of reaction at minimum temperature is obtained. Results show that platinum is a perfect catalyst for converting H₂S to H₂ and it has a perfect performance (98%) after 5 μs at low temperature of 227°C.     

Exhaustive search of the configurational space of heat-shock protein 90 with its inhibitor by multicanonical molecular dynamics based dynamic docking

Multicanonical molecular dynamics based dynamic docking was used to exhaustively search the configurational space of an inhibitor binding to the N-terminal domain of heat-shock protein 90 (Hsp90). The obtained structures at 300 K cover a wide structural ensemble, with the top two clusters ranked by their free energy coinciding with the native binding site. The representative structure of the most stable cluster reproduced the experimental binding configuration, but an interesting conformational change in Hsp90 could be observed. The combined effects of solvation and ligand binding shift the equilibrium from a preferred loop-in conformation in the unbound state to an α-helical one in the bound state for the flexible lid region of Hsp90. Thus, our dynamic docking method is effective at predicting the native binding site while exhaustively sampling a wide configurational space, modulating the protein structure upon binding.